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ABSTRACT 


SYSGEN is a production costing and reliability model of electric utility 
systems. Hydro-electric, storage, and time-dependent generating units can be 
modeled in addition to conventional generating plants. 

This user documentation describes the SYSGEN model and its links with 
other JPL simulations. Input variables, modeling options, output variables, 
and report formats are explained in detail and illustrated with examples. The 
appendixes contain sample computer output and instructions for running the 
program on an IBM batch system. SYSGEN also can be run interactively by using 
a program developed at JPL called PEPS (Front End Program for SYSGEN). A 
format for SYSGEN input variables which is designed for use with FEPS is 
presented. 


FOREWORD 


SYS6EN is a computer program which simulates the production costs and 
reliability of an electric utility system with and without time-dependent 
generating units. It is one of several simulations used at JPL for 
calculating the breakeven costs of photovoltaic and solar thermal electric 
systems in grid-connected applications. 

JPL received SYSGEN in 1979 from the Massachusetts Institute of 
Technology (MIT). Significant modifications have been made in the SYSGEN 
program since then, based on the work of Dr. Donald Ebbeler and Dr. George 
Fox, both of JPL. These modifications, listed in Appendix C, do not affect 
the inputs to the model; thus, this user documentation may be used with both 
the MlT and JPL versions of SYSGEN. 
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SECTION I 


INTRODUCTION 

A. THE USE OF COMPUTER SIMULATIONS IN THE BREAKEVEN COST CALCULATION 
SYS6EN is one of several simulations used in the breakeven cost 
calculation for photovoltaic and solar thermal electric systems. Together, 
these simulations incorporate sufficient technological and financial detail to 
allow realistic estimation of the value and effects of solar generation to a 
grid system. The rest of this section explains how the JPL simulation models 
calculate these values. 

Four major areas are modeled: 

1) the lifetime cost and performance of the solar plant 

2) the production cost and reliability of electric power generation by 
the utility, with and without the solar plant 

3) optimal capacity expansion planning for the utility over time, with 
and without the solar plant 

4) the financial structure and environment of the utility. 

Electric power systems are assumed to operate to meet a fluctuating 

power demand at minimum cost.* The demand for electricity generally follows a 
predictable daily pattern, with peak periods occurring around noon and early 
evernng. Utilities adjust their generating plant mix according to these 
changes in demand. Generating plants with low operating costs are run all the 
time to meet the minimum power demand and are called base loaded units. 

Utilities in reality are constrained by electric stability requirements 
imposed by the transmission network. However, a complete model of 
operating a power system would require detailed models and data for each 
generator and transmission Tine and would be far too complex for our 
purposes. Most production costing models, including the one used by JPL, 
do not consider transmission or stability constraints. 


Plants with high operating costs but low capital costs are brought on line 
during the peak hours and are referred to as peak loaded units. Thus» as a 
utility Increases its output, the marginal cost of generation also increases 
as more and more expensive units are brought on line to meet demand. 

The value of a solar plant to a grid Is defined in terms of fuel 
displacement and capacity displacement. Displaced fuel is fuel saved because 
solar generation is used to meet load instead of using an existing 
conventional plant. Capacity displacement occurs when the addition of new 
conventional capacity can be deferred or cancelled if a solar plant is added 
to the grid. 

A utility can reduce its operating costs either by reducing the size of 
the peak loads (e.g., with time-of-day pricing) or by displacing expensive 
peak unit generation with cheaper limited energy sources, such as solar, 
storage, or hydro plants. It should be apparent that the value of a solar 
plant will depend upon what kinds of generating units it displaces. The value 
of the fuel displaced depends upon the correlation during the day between 
solar output and the load on the utility: If the bulk of solar output occurs 

during off-peak periods, then cheap base loaded plants are displaced rather 
than expensive peaking units. It is for this reason that utility and solar 
plant cost and performance simulations are done on an hourly basis rather than 
averaged over days or months. Averaged data will not correctly capture the 
difference between the value of peak and off-peak generation and thus will 
underestimate or overestimate the displaced fuel value attributed to the solar 
plant. 

The lifetime cost and performance (LCP) simulation of a photovoltaic 
plant performs four functions: 1) it describes initial design and 
construction of the power plant, 2) simulates hourly photovoltaic system 
power output, 3) analyzes the long-term effects of exposure to an outdoor 


environment and operations/maintenance policies, and 4) finds the capital and 
recurring production costs of operating the system. The LCP model contains 
both deterministic and probabilistic characteristics! Module failure is 
modeled using transition probabilities; explicit time-dependent features are 
modeled deterministically,* 

The cost and performance of solar thermal systems are estimated by the 
Solar Energy Simulation model (SES). The SES model comprises three programs 
which perform the following functions; 1) evaluation of collector field 
perrormance for specific insolation and temperatures, 2) calculation of the 
performance of a fixed-rated power plant with different collector and storage 
sizes, and 3) calculation of the minimum energy cost of the plant. These 
models are coupled with the SYSGEN utility grid simulation to determine the 
effects of the solar plant upon grid operating cost and reliability. 

The value of the solar plant is also determined by the timing of its 
introduction to the grid system. The optimal time for bringing solar capacity 
on-line will be determined in part by future availability of fuels and costs 
of conventional power plants. JPL's simulation models allow the testing of 
different assumptions about the above factors. The production costing model 
(SYSGEN) can be run for up to a 35-year time period, incorporating different 
capacity expansion plans, cost escalation rates, and changes in customer 
demand. The results of these different production cost scenarios are then 
combined with investor specific financial and Investment data to determine how 
the feasibility of investment in the solar plant changes over time. 

Utilities plan their generation mix to meet a certain reliability 
requirement. One measure of reliability is called the loss of load probability 

* C.S. Borden, ’'Computer Simulation of the Operations of Utility Grid 
Connected Photovoltaic Power Plants," paper presented at The National 
Computer Conference, May 19, 1980. 
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(LOLP), which may be defined as the probability that the utility cannot 
generate enough power to meet demand at any point in time, either because of a 
sudden surge in demand or an unexpected plant failure. On an hourly basis, 
the LOLP allowed is usually specified as a failure to meet load for one hour 
out of every ten years.* In order to meet this reliability criterion, 
utilities keep additional "unnecessary" units on spinning reserve, units which 
are running and can be connected to the grid immediately if needed. The 
production costing model uses a probabilistic method to incorporate 
uncertainty due to conventional plant failures and calculates the LOLP of the 
system. Because of this LOLP constraint, it is necessary to know the effect 
of a grid-connected solar plant on the reliability of the system. Solar 
plants differ from conventional units in that uncertainty in their performance 
is due not only to mechanical failure but also to the vagaries of the 
weather. Unfortunately, if scenarios are modeled using historical weather 
data (as they are at JPL), neither outages due to weather nor due to 
mechanical failure are included in the LOLP calculation for the grid system. 
This does not affect the estimates of the fuel displacement value for solar, 
but it presents some problems in planning capacity displacement. 

One approach to this problem is to model both utility load and solar 
output as random variables, whose distributions are dependent on other 
measurable variables such as temperature, insolation, and time. This method 
would work particularly well in southern California, for example, where summer 
peak demano for electricity is highly dependent on the temperature because of 
the demand for air conditioning. 

More specifically, in a period of ten years, the expected number of 

hours for which there is a capacity deficiency to meet the hourly load 

demand is one. 
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The capacity displacement attributed to a solar plant is estimated with 
the optimal capacity expansion and the production costing/reliability models. 
The capacity expansion program (SCYLLA) uses a linear programming algorithm to 
minimize the cost of the generation mix for given changes in demand and 
production costs. Optimal changes in generation mix over time with and 
without solar plants are determined and are then tested for their 
reltabillty. If the addition of conventional capacity can be deferred when 
solar plants are added, then a capacity credit for the deferred capital costs 
is attributed to the solar plant. 

The financial analysis model sfAPSEAM) incorporates 1) the estimated 
capital and production costs of both the solar plant and the utility, 2) the 
benefit stream to the utility in terms of fuel and capacity displacement over 
the lifetime of the solar plant, and 3) the optimal capacity expansion plan 
for the utility. These cost and benefit streams then can be tested in various 
financial environments, which include consideraton of inflation, escalation 
and interest rates, depreciation, and taxes to determine the net present value 
of an investment in the solar plant. The breakeven cost of the systein is that 
cost which gives a net present value of zero for the investment. 

Electric utilities use similar computer models for determining daily 
dispatch strategies and capacity expansion plans. Thus, our modeling 
procedures should be accepted as valid by the utility industry. The inter- 
faces of the simulation programs are shown in Figure 1. 
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B. THE SYSGEN USER PACKAGE 

The SYSGEN computer simulation is a production costing and reliability 
model of an electric power system. It may be used to study the effects of 
different assumptions about customer demand, fuel costs, or utility 
characteristics on system cost and reliability. 

This package is designed for users who are not conversant with computer 
programming or utility engineering. Thus, descriptions of the actual modeling 
procedures are brief and necessarily general, and may not be detailed enough 
for some users. A co)nplete list of the technical documentation for SYSGEN may 
be found in the references of this paper. 

The model consists of three computer programs. The first, named 
ELECTRA, calculates a load duration curve (LDC) from hourly load curves, and 
models time-dependent power plants as increases or decreases in the net load 
on the system. Two LDCs are produced: one with the time-dependent generation 

and one without. The second program, SYSGEN, uses these LDCs and data 
describing the power plants in the system to calculate system cost and 
reliability with and without the time-dependent units. SCYU.A, the third 
program, uses these results to calculate the fuel and capacity credits of 
the time-dependent generation to the power system. The diagram in Figure 2 
shows the program interactions. 

Sections II through IV of this report explain the inputs, outputs, and 
modeling procedures used in the three computer programs. Section V contains 
worksheets which may be used to set up the inputs to the programs. These 
worksheets correspond directly to an interactive computer program named PEPS 
(Front End Program for SYSGEN). FEPS may be used to create computerized data 
bases for these simulations, and is explained below. 
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Figure 2. SYSGEN Model Overview 




C. FEPS— THE FRONT END PROGRAM FOR SYS6EN 

FEPS is an Interactive computer program which may be used to create and 
modify data bases for SYSGEN simulations. The program transforms inputs into 
the proper format and files and adds the JCL (Job Control Language) necessary 
for SYSGEN, ELECTRA, and SCYLLA runs. It is presently set up on the 
California Institute of Technology (CIT) VAX time-sharing system, and may be 
accessed either from terminals at CIT or from remote terminals at JPL. 

FEPS was written to be used in conjunction with this user package. 
Documentation for FEPS is available upon request. Data bases for SYSGEN still 
may be created from fixed format punched cards if so desired, and the OCL for 
running SYSGEN with card inputs is explained in Appendix B. The input formats 
for cards are fully described in the original MIT documentation for SYSGEN, 
also available upon request. 


SECTION II 
ELECTRA 


A. THE LOAD CURVES 

ELECTRA models time-dependent plants as decreases in the net load on the 
utility grid system. Up to four hourly load change curves may be input in 
addition to the original custcmer demand curve. The hourly load curves are 
converted into load duration curves that represent the net load with and 
without the time-dependent generation. 

There may be up to 8,784 values in the customer demand curve and in each 
of the load modification curves. These values do not have to be hourly, but 
the values for demand and modification must match, i.e., they must start at 
the same point in time, the time interval used must be the same, there must be 
the same number of points in each curve, and the energy units must be the same 
(MW or kW). 

The load duration curve is determined by the percent of time that load 
is greater than or equal to a percentage of peak load. Specifically, the 
dependent variable (Y axis) is the 

percent of time demand > X • PEAK LOAD/n, 

where X = 1, 2, 3, n, and n - number of points in the load shape. 

The "number of points in load shape" chosen determines the size of the 
MW intervals used in approximating the curve: The smaller the interval size, 

the more closely the true load duration curve is approximated. The user must 
specify the peak of each original load curve and the number of points in the 
load shape to be used in calculating the load duration curves. An example 
should be enlightening. 
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Suppose peak demand is 60 MW. We want to km ^ the percent of time that 
the load is greater than or equal to Y MW, 0 < Y <60, where Y = X • 60/n, 
with n and X defined as before. Table 1 illustrates the calculation of the 
load duration curve using three different values of n. (See also Figure 3.) 

Table 1. Determining the Load Duration Curve 



Number of Points in Load Shape 


5 Points 

10 Points 

20 Points 


of Time 

Load is >1 X 60/5 = 12 

2 X 60/5 = 24 

3 X 60/5 = 36 

4 X 60/5 = 48 

5 X 60/5 = 60 


% of Time 

Load is ^ 1 X 60/10 = 6 

2 X 60/10 s 12 

3 X 60/10 = 18 

• . 
t 

10 X 60/10 = 60 


% of Time 

Load is ^ 1 X 60/20 = 3 

2 X 60/20 = 6 

3 X 60/20 = 9 

♦ 

i 20 X 60/20 = 60 MW 


Thus, if there are 5 load shape points, the load interval will be 12 MW; 10 
load shape points gives a 6 MW interval, and so on. It should be noted that 
the cost of ELECTRA runs increases almost exponentially as the number of 
points in the load shape goes above 40. 

Frequency curves are also formed from the original and modified load 
curves. The frequency curve is determined by the number of times the load 
rises across a given MW level in time period T, normalized by the number of 
hours in T. The given MW levels are the same as those chosen for the load 
duration curve. 



Up to 52 frequency and load duration curves may be formed from the 
original load curve, e.g., given a year's worth of hourly load data, the user 
can model each week of that year separately (or each month or season). The 
number of curves desired and the number of hours in each must be specified by 
the user. 

The reliability of the utility system Is estimated by adjustinq the load 
duration curve to account for possible plant failures. The adjusted LDC's are 
called equivalent load duration curves. These curves are actually calculated 
In SYS6EN, but their description here Is a logical extension of the discussion 
of the load duration curve. 

If a generating unit may be loaded Incrementally (1.e., if it has valve 
points), then the unit may be treated as either a single Increment or as 
multiple increments In the calculation of the equivalent load duration curve. 
If the former Is chosen, then the unit Is either all on or all off; If the 
latter is chosen, then partial outages may be modeled. 

Figure 4 inustrates the simpler case of treating the unit as a single 
increment. Suppose two plants are being loaded. There is some probability 
that the first plant will fail and the second plant will be faced with the 
entire original load curve. If the first plant does not fail, then the second 
plant will be faced with the original load minus the capacity of the first 
plant. The equivalent load duration curve when the second plant Is loaded Is 
the sum of the curves weighted by their respective probabilities. 

The equivalent LDC In Figure 4 was derived as follows. The first plant 
to be loaded has a capacity of 5 MW. Let be the state in which the 5 MW 
unit does not fail, and $2 be the state In which the 5 MW unit does fail. 

The probability of S-j occurring is 0,9; the probability of $2 occurring is 
0.1. From the figure it can be seen that a load level of 20 MW occurs 0% of 
the time in S|, but in $2 it occurs 501^ of the time. Then the weighted 
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probability that a load of 20 MW occurs Is: 0.9 (0<) + 0.1 (5055) * 55^. This 
1s point C on the graph. A load of 10 MW occurs 1005J of the time In both 
and Sg, so the weighted probability Is 0.9 (100%) *♦ 0.1 (100%) » 100%, 
giving us point A. The weighted probabilities of 17 and 22 MW loads were 
derived In the same fashion: For point B, we have 0.9 (30%) + 0.1 (80%) * 
35%, and for point D It Is 0.9 (0%) + 0.1 (30%) * 3%. 


B. TIME PARAMETERS 


SYSGEN’s time parametefs may be set to model many different situations. 
There are three basic units to specify: the number of periodsi subperiods In 
a period, and weeks in a subperiod. A period Is usually Interpreted as a 
year, and a subperiod may represent a week, a month, or any fraction of a 
year. The time units to be specified are as follows; 


NUMBER OF PERIODS 

"HOURS'* IN A TIME PERIOD 


SUBPERIODS 


1 ) 


"WEEKS" 


"HOURS" IN A "WEEK" 


An Integer value; 1 < P < 34. 

Each time period will have the 
same number of "hours"; this need 
not be an integer value; 

1 < H < 8,784. 

An integer value; 1 :< SP S 52, 

Each period will have the same 
number of subperiods. 

An Integer value; 1 < W < 99. The 
number of "weeks" per subperiod 
may varied . 

Each "week" will have the same 
number of "hours"; this need not 
be an Integer value. 


To illustrate this, consider the following three examples: 


a) 

b 

c) 

d 


e) 

2) a) 
b) 

5i 

e) 


number of periods * 10 (10 years) 

"hours" in a period = 8736 (24 hours x 364 days) 

number of subperiods = 4 (seasons) 
number of "weeks" In subperiod 1 ~ 5 
number of "weeks" in subperiod 2 = 2 
number of "weeks" in subperiod 3^3 
number of "weeks" in subperiod 4 * 2 
"hours" per "week" = 720 (24 hours x 


(winter months) 
(sprino months) 
(summer months) 
(fall months) 

30 days) 


NOTE THAT THE 
"WEEK" VARIABLE 
EQUALS A MONTH 
IN EXAMPLE 1 


number of periods = 5 (5 years) 

"hours" in a period = 2,920 (data is measured at 3-hour Intervals 

for 365 days) 

number of subperTods> 2 (winter and summer) 

number of "weeks" in subperiod 1 - 20 (20 "weeks" in winter) 

number of "weeks" in subperiod 2 = 32 (32 summer weeks) 

"hours" in a "week" = 56 (data measured at 3-hour intervals 

for one week) 
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3) a) number j)f periods « 12 (12 months) 

b) "hours'* in a period * 672 (hours per 4-week month) 

cj number of subperiods * 4 (4 weeks per month) 

dj "weeks" In each subperiod * 7 (7 days per week) 
e) "hours" per "week" »• 24 (24 hours per day) 

These examples should show that practically any time frame may be used, 
as long as the variable values are internally consistent. Example 2 
illustrates that the "hour" variable must correspond to the time value of the 
original load datai In this instance, one "hour" equals three hours. In 
examples 1 and 3, the "week" variable is equivalent to a month and a day, 
respectively. Please note that the assigned values of the "hour" and "week" 
variables will be used throughout the program. For example, a later input to 
the program is "the number of weeks per subperiod a generating plant is down 
for maintenance." Then if the "week" variable is one day and the unit will be 
down for 7 days, the "number of weeks per subperiod for maintenance" will be 
7. Other variables which are input in "weekly" values are the weekly energy 
capacity of storage and hydro units and the installment and retirement weeks 
of generating units. 


SECTION III 
SYS6EN 


The objective of this program Is to find an operating schedule which 
minimizes production costs, and to estimate the frequency, duration, and 
probability of loss of load for a given mix of generation units and a given 
customer demand. The program uses a modified Booth-Baler laux technique which 
t’eats plant outages as randomly occurring loads on other plants In the 
utility system.* 

The Input data required for SYS6EN are divided Into the following 
sections: 

A. Economic parameters 

B. Program operation options (Input options) 

C. Transmission and distribution losses 

D. Generating class data 

E. Generating unit data 

F. Loading order 

G. Report optit ", and outputs from SYSGEN 


This methodology is described in S. Finger, "Electric Power System 
Production Costing and Reliability Analysis Including Hydro-electric, 
Storage, and Time Dependent Power Plants," MIT Energy Lab Technical 
Report ImIT-EL- 79-006, February 1979, and in D. H. Ebbeler, "Electric 
Power Generation Systems Probabilistic Production Costing and 
Reliability Analysis," (JPL working paper), April 1981. 
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A. ECONOMIC PARAMETERS 

The production costs in each time period may be reported either as 
nominal or present values. The economic parameters used by SYS6EN are 1) the 
overall inflation rate, 2) the real discount rate, and 3) the nominal 
escalation rates for fuel and operations and maintenance (OiM) costs.* These 
values are entered as fractions, so a lOilS discount rate is input as ”0.10.” 
The inflation and discount rates are assumed to be constant throughout the 
simulation. The escalation rates, however, may vary by time period. All 
rates are assumed to correspond to the designated time period length and are 
constant over the subperiods. 

Escalation rates may be varied by generation class type, but not by 
individual plants. These input data are thus entered with the generating 
class data rather than with the other economic parameter data. 

All input costs, such as fuel and cost per start-up, should be inpu.t in 
the same year's nominal dollars. The user may choose in what year's dollars 
the production costs will be reported. If present worth values are not 
desired, the user must set the discount rate to zero. 

It Is important to note that if the year the input costs are in is 
before the year the simulation is started, then all costs will be adjusted to 
start year nominal dollars using the general inflation rate, not the 
product-specific escalation rate. The easiest way to avoid confusion here is 
to input the costs in the start year's nominal dollars. 

The following formulas are used to calculate the present worth factor 
and the esvralation factors for fuel and O&M costs. 


An escalation rate is basically a product-specific inflation rate, e.g., 
a rise in fuel costs alone Is designated by an escalation rate which, 
these days, is higher than the general inflation rate. 


Let 

T 

e 

current time period 



I 

sc 

generating class 



DR 


discount rate 



CPI 


consumer price index 



ERom(I.T) 

n 

escalation rate for O&M costs for class I 

in time period T 


ERf(I,T) 

sc 

escalation rate for fuel costs for class 
T 

T in time period 


GONVRT 


conversion factor from input year nominal 
report year nominal dollars 

dollars to 


CiN 


cost in the input year nominal dollars 



PWFx 

£ 

present worth factor in time period T 



£ 

(1 nR)-T 


then 





ESCOM 

j (I) 


escalation factor for n?iM for class 1 at 

time period T 



£ 



ESCFx (I) 

£ 

escalation factor for fuel for class I at 

time period T 


TT [l + ERp (i,t)] 

t=i ' 

and 


Cj = present value in report year dollars of a cost incurred in 

time period T of the study 

Q . g M BI . eSCt (I) / PWFr 

n+r.pi)‘ ' 

( ESCOMt (I) for an O&M cost 

where ESC (I) - { ' 

( ESCFj(I) for a fuel cost 
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B. INPUT OPTIONS 

SYS6EN has a set of logical variables that can be used to control the 
level of detail in the simulation. These options will supersede any input 
parameters. For example, if the spinning reserve option is set to false and 
the spinning reserve requirement is set to 200 MW, no spinning reserve will be 
used in the simulation. 

Not all of SYSGEN's input options are explained in this section. The 
options on hydro and storage dispatch, time-dependent plants, and loading 
order are explained in later sections as appropriate. 

1. Spinning Reserve Option; MSP IN 

The spinning reserve algorithm is implemented after the economic loading order 
is set up. The loading order is modified, if necessary, to meet the input 
spinning reserve requirement. Three variables are needed to specify spinning 
reserve: 

1) The required reserve, stated in one of the following three ways: 

a) as a percent of peak load 

b) as an absolute MW value 

c) as a fraction of the largest unit on line. 

2) The maximum percent of any unit to be credited to spinning reserve. 

3) The maximum number of units to be displaced in the economic loading 

order to meet spinning reserve. 

In reality, no utility ever loads its plants as SYS6EN does, i.e., starting 
cold with no energy generated and having to bring base units on-line. In a 
real utility, the base units are always running and the amount of capacity on 
spinning reserve is relatively constant. But as a computer model, SYS6EN does 
not see a unit as available for spinning reserve until it is partially loaded. 
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Thus, SYS6EN calculates the available reserves as the remaining capacity of 
units that have been partially loaded, up to the maximum percent of the unit 
specified in (2) above.* The following scenario should serve to illustrate 
this. 


Suppose the loading will be started with two coal plants having 
capacities of 400 and 600 MW. Each coal plant has five loading increments, as 
follows. 


Incremental 
Unit Capacity 

Coal 400 100 

60 
80 
80 
80 


Incremental 
Unit Capacity 

Coal 600 150 

90 
120 
120 
120 


The spinning reserve has been set to 300 MW. A maximum of 50J5 of any 
unit m^y be credited to spinning reserve. Further, assume that the economic 
loading order is to load all of the 600 MW unit and then the 400 MW unit. It 
can be seen, however, that if all of the 600 MW plant is loaded before the 400 
MW plant, the available reserve will be less than 300 MW for some period of 
time. SYSGEN therefore will load the first two increments of the 600 MW unit, 
then load the first increment of the 400 MW unit. The "load/not load" 
decisions made by SYSGEN for the first few increments are shown in the table 
below. Note that the available reserve is not allowed to fall below 300 MW 
and that no more than 50% of a unit is ever credited to spinning reserve. 


This method of calculating the available spinning reserves is valid in 
the deterministic case, but it may not be appropriate when plants are 
modeled probabilistically since a zero forced outage rate is implicitly 
assumed for the plant when it is on spinning reserve. 
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Table 2. The Effect of a Spinning Reserve Requirement 


on Loading Or^er 


Unit 

Increment 

Number 

Increment 

Capacity 

Reserves 

Available 

System 

Capacity 

Load 

Oecision 

600 

1 

150 

300 

150 

Yes 

600 

2 

90 

300 

240 

Yes 

600 

3 

120 

240 


No 

400 

1 

100 

500 

340 

Yes 

600 

3 

120 

440 

460 

Yes 

600 

4 

120 

320 

580 

Yes 

600 

5 

120 

200 


No 

400 

2 

60 

320 

640 

Yes 


The loading order option (UORDOP, explained on p. 39) will supersede the 
spinning reserve requirement. For example» suppose the loading order option 


chosen required that all base units are loaded before any intermediate, and 
all intermediate units are loaded before any peak units. Then no intermediate 
units can be brought up to satisfy the spinning reserve requirement if not all 
base units have been loaded, no matter what value of spinning reserve is input. 

2. Multiple Increment Option; MULT 

MULT controls the multiple increment loading algorithms. If MULT is set 
to false, units are modeled as on/off variables (i.e., as a single 
increment). The single increment characteristics can be input or, if left 
blank, will be computed from the data for multiple increments. Note that if 
MULT is set to false, then MSP IN is automatically set to false by the program 
because of the method used to calculate spinning reserve. 


3. Frequency Curves Option: MFREQ 

MFREQ controls the frequency and duration algorithms. If MFREQ is set 
to false, no frequency curves are read in, and no frequency calculations, such 
as the expected number of startups, are made. 

SYSGEN uses the average duration of each load level in calculating the 
amount of capacity available for spinning reserve. If MFREQ is set to false, 
but the spinning reserve option, MSPIN, is set to true, then the frequency of 
every load level will automatically be set to 1/number of "hours” in a day. 
This is a modification of the original MIT SYSGEN program. If MFREQ is false 
in the MIT version, then the frequency of every load level is assumed to be 
one. 
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c. 


TRANSMISSION AND DISTRIBUTION (T&D) LOSSES 

T&D losses may be simulated for both conventional and time-dependent 


plants. Conventional plant losses are accounted for by dividing the energy 
output by a penalty factor specified by the user. Thus, a penalty factor of 
1.0 will mean there are no T&O losses for that plant. (Note that the penalty 
factor cannot be less than 1.0.) 

Time-dependent plant losses are modeled as decreases in the load 
reduction curves. The marginal losses in MW at up to ten load levels may be 
specified. The load levels are given as percentages of peak demand, and the 
intervals between demand levels do not have to be equal. Linear interpolation 
is used to find the intermediate losses. 



D. GENERATING CLASS DATA 

Each plant in SYSGEN must be assiqned to a qenerating class, such as a 
hydro, diesel, or storage class. The Input variables which are specified by 
class type are the cost escalation rates, the forced outage multipliers 
(explained below), and the generation type, i.e., base, intermediate, or 
peaking generation. Up to 34 class types may he created. 

Hydro, storage, and time-dependent plants must be assigned to 
pre-spec if ied SYSGEN class types: 

1) All plants must be assiqned to class type *'HYD0'* when using 

Ftps, and to "CHY" when card input is used. 

2) All storage plants must be assiqned to class type "STOR" when using 
FEPS, and to "STO” when using cards. 

3) All time-dependenH plants must be assigned to class type "LORD" 
(short for "load reduction") when using FEPS, and to "TOP" when 
using cards. 

If these plants are not assigned to the proper classes, then SYSGEN will treat 
them as conventional units. 

The forced outage rate of a plant is the percent of time the plant 
breaks down and is not operating (sometimes called "unscheduled maintenance"). 
A forced outage rate is input for each plant in the simulation, A variable 
called the "forced outage multiplier" may be used to increase a plant's forced 
outage rate during specific time periods. This multiplier would be used, for 
example, to account for variations in the reliability of a newly installed 
plant. 

The forced outage multiplier is used as follows: 

Let FOR-j s input forced outage rate for unit j. 

IFOR-j-j = forced outage multipler for unit j in 
year i. 
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then 


FORj X IFOR^j » forced outage 
rate for unit j in year i used in 
SYS6EN calculations of LOLP. 


Two things should be noted here. The first is that forced outage 
multipliers are assigned by generating class, not by plant. Thus, if you only 
want to vary the failure of one plant, you must create a new class for it. 

The second point is that SYSGEN includes the time a plant is out on scheduled 
maintenance in the calculation of the plant's effective capacity.* Thus, 
forced outage multipliers should not be used to account for time that the 
plant is down on scheduled maintenance. 

The escalation rates and forced outage multipliers are entered in the 
form of "tables" (vectors, actually) with rates assigned by class type for one 
or more periods in the simulation. A maximum of 10 escalation and 10 outage 
tables for 10 periods (values) each may be specified. An example may he of 
use. 

Suppose a simulation runs for three periods and has two plant classes; 
steam turbine and gas turbine. If the fuel and O&M costs for these classes 
all escalate at different rates, four escalation tables must be input. These 
might be 


Glass Type 

One 

Period 

Two 

Three 

Table Number 

Steam Turbine 

Fuel 

0.06 

0.07 

0.08 

1 

O&M 

0.05 

0.05 

0.06 

2 

Gas Turbine 

Fuel 

0.08 

0.10 

0.11 

3 

O&M 

0.06 

0.07 

0.07 

4 


The effective capacity of a plant is defined as the capacity that a 
100?t-re liable plant would have to have in order to generate an equal 
amount of power over the time period. 
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On the other hand, O&M costs for both classes mioht increase at the same 


rate. Then only three tables need be specified: 


Class Type 

One 

Period 

Two 

Three 

Table Number 

Steam Turbine 
Fuel 

0.06 

0.07 

0.08 

1 

Gas Tut bine 
Fuel 

0.08 

0.10 

0.11 

2 

O&M 

0.05 

0.06 

0.04 

3 


If the number of escalation values entered is less than the number of 
time periods in the study, the last value entered is assumed to hold true for 
the remaining time periods. Thus, if an escalation rate will he constant in 
the study, that value need only be entered once in the first time period. 

Note that the same escalation rate table may be assigned to any number 
of classes and for either fuel or O&M costs. 

The forced outage multipliers are set up in the same way, with one 
exception. If a value in the multiplier table for a given period is set to 
zero or left blank, the forced outage rate for plants in that class and time 
period will be calculated as zero (as can be seen from the eouation on 
page 28: ^ ^ forced 

outage rate will not be altered for a class, a value of "1.0" should be 
entered in the multiplier table. Once a value of 1.0 has been found in that 
table, all subsequent periods in that table are assumed to have the value 1.0 
unless other values are explicitly entered. Consider the following three 


multiplier 

tables: 

Period 



Table 

1 

2 

3 

4 

1 

1.2 

1.1 



2 

1.2 

1.1 

1.0 


3 

1.0 

i* 

• 
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With Table 1, the forced outage rates caTculated will be higher than 
originally set in time periods 1 and 2, but in time periods 3 and 4 the rates 
will be zero (thus simulating a totally reliable unit). Table 2 will cause 
the program to apply the original forced outage rate in time periods 3 and 4, 
with higher rates in the first two periods. If Table 3 is assigned to a 
class, no forced outage multiplier will be used, so the original forced outage 
rate is applied in all four time periods. 

In summary, the generation class data required is 

1) Number of generation classes in the mix (e.g., 2 nuclear plants, 4 
coal fired, and 3 gas turbines = 3 generating classes) 

2) Class name of each (e.g., COAL, HYDO) 

3) Whether the class is base, intermediate, or peaking generation 

4) Cost escalation tables for fuel and OiM 

5) Forced outage multiplier tables 


E. GENERATING UNITS 

There are three special types of generating units in SYS6EN: hydro, 
storage, and time-dependent plants. All other generating plants are called 
conventional units. 

The following list is a sumnary of the input data required for all types 
of generating units, Hydro, storage, and time -depen dent plants require 
additional inputs which are explained at the end of this section, 

1) The generation class to which the unit belongs 

2) Incremental loading data , if relevant 

3) Installment year and week (e.g., a unit installed on February 2, 
1979, would be year » 1979, week = 5, if real weeks are used) 
Retirement year and week (e.g., a unit retired on December 1, 1980, 
would be year « 1980, week » 48, as above) 

Fuel cost for the unit ($/MBtu) 

6) Variable O&M cost ($/MWh) 

7) Cost per startup ( S/startup) 

8) Cost per MWh to keep unit as spinning reserve without generating 
power ($/MWh) 

9) Mean time to repair after failure (hours) 

^0) Penalty factor , used to account for unit-specific transmission 
losses 

11) Preventive maintenance data for each unit. It is important to note 
that the maintenance cycle is referenced to the installation date 
of the unit, not the start of the simulation. 
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It is assumed that maintenance is cyclically scheduled, and the 
cycle may be from one to ten periods long. Within each cycle, the 
subneriods during which maintenance is performed and the number of 
weeks in those subperiods the unit is unavailable due to 
maintenance must be specified. The number of weeks is an integer 
value which may vary by subperiod. If a plant is out more than one 
subperiod in a row, this must be entered explicitly. 

Example: Simulation is started on January 1, 1988, and the unit is 
installed January 1, 1989. It is being put on a two-year 
maintenance cycle, with the following schedule: 


June 20, 

1989 

1 week 

December 

10, 1989 

2 weeks 

June 15, 

1990 

2 weeks 

December 

n, 1990 

3 weeks 


Period length = 1 year, 12 subperiods (months) per period. Then 


cycle length = 2 and 

Subperiod Weeks 

Period 1: 6 1 

12 2 

Period 2: 6 2 

12 3 
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Capacity , heat rate (MBtu/MWh), and forced outage rate (average for 
the plant and for each loading increment of the plant, if 
applicable) 

For example, a 400 MW plant with four loading increments might have 
the following characteristics; 


Increment 

1 

2 

3 

4 


Cumulative 

Capacity 

80 

160 

280 

400 


(MBtu/MWh) 

16 » 9 
12.2 
11.8 
12.0 


Forced 
Outage Rate 

0.03 

0.0194 

0.019 

0.0279 


The forced outage rate is actually a conditional probability. Each 
increment has an independent probability of failing, but increments 
must be loaded sequentially. Thus, the forced outage rate of 
increment i is the conditional ^irobability that all increments 
before i have been loaded and increment i cannot be loaded. In the 
example above, each increment had the following independent 
probabilities of being loaded: 


Increment 

1 

2 

3 

4 


Probability 
of Failure 

0.03 

0.02 

0.02 

0.03 


Probability 
of Success 

0.97 

0.98 

0.98 

0.97 


The probability that cumulative capacity equals zero is simply the 
probability that the first increment fails, so the forced outage 
rate for the first increment is 0.03. In order for cumulative 
capacity to be exactly 80 MW, the first increment must be loaded 
and the second increment must fail. Then the forced outage rate 
for the second increment is (0,97)(0.02) = 0.0194, and so on. 

In practice, the incremental probabilities are not available. The 
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forced outage rates given are usually a total outage rate for the 
plant. In the example above, the total forced outage rate is 
0,0963 (i.e,, the plant will be be generating at full capacity 
90.37< of the time). This total outage rate may be obtained two 
ways: 1) It is the sum of the incremental forced outage rates, or 
2) it is one minus the product of the independent success 
probabilities. 

1. Hydro Units 

The only additional data required for a hydro unit is the "weekly” 
size of the hydro reservoir in each subperiod, expressed in megawatt hours. 

Hydro power can be dispatched in two different ways. The first is to 
delay loading until the hydro unit can generate all its energy at full 
capacity. This delay is determined by the equivalent demand curve: The hydro 

unit will not be loaded until the remaining demand for energy is less than the 
available hydro energy. The second dispatch strategy is to simply load all 
hydro units first at a reduced capacity to generate all of their energy. If 
the input option MOLAY is set to true, then the first hydro dispatch is used. 
If MDLAY is set to false, then the hydro units will be loaded first as base 
generation, even if the hydro units are classed as peak or intermediate 
generators. 

The first strategy is expected to give a more cost efficient use of the 
hydro energy than the second strategy, since more energy is removed from the 
peak when the hydro dispatch is delayed.* 

A second dispatch option, MOVE, allows SYS6EN to interrupt the loading 

~* The first strategy is optimal in deterministic production costing, but 
may not be optimal for probabilistic production costing. See D. H. 
Ebbeler, "Electric Power (Generation Systems Probabilistic Production 
Costing and Reliability Analysis," (JPL working paper), April 1981, and 
G. Fox, "A Stochastic Formulation of Electric Generation System 
Reliability Measures," (JPL working paper), January 1981, 
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of a conventional unit and back It off until the total energy demand exactly 
equals the hydro energy available. If MOVE is not true, then hydro units may 
only be loaded after increments of conventional plants have been fully loaded. 

Purchased power can be. modeled in SYSGEN as a hydro plant. The only 
difference is that the fuel cost in $/MBtu should be set to the purchase price 
in $/MWh, and the "heat rate" of the plant is set to 1,0. 

2, Storage Units 

Two distinct steps are involved in modeling storage units; 
charging the units with excess energy and then dispatching the stored energy. 
The following data are required for each storage unit: 

1) The charging capacity, in megawatts 

2 ) The forced outage rate for the charging cycle 

3) The overall efficiency of the storage and generation process, 
measured by the product of the charging and discharging 
efficiencies (or equivalently, energy available from storage/energy 
generated to charge storage) 

4) The reservoir size in megawatt hours available from storage in a 
"week." This value should be the most energy available from 
storage to meet demand when charging and discharginq efficiencies 
are taken into account. 

The input option MSTOR determines how storage is handled. If MSTOR is 
set to true, then the amount of energy stored is determined by the expected 
excess energy available from base units, expected storage plant failures, and 
the expected cost of the stored energy. If there is more than one storage 
unit, then they are ranked so that the largest ones are filled first. It 
should be pointed out that these storage algorithms are not optimal. The 
program assumes that the storage units will be filled, and then asks when they 
should be loaded, missing the question of how much energy should actually 
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be stored, based on the marginal cost of adding energy to storage. 

If MSTOR is set to false, then the energy available from storage is 
taken as the input storage size. The marginal cost of storage is similarly 
set to the average cost of base loaded energy. 

The input options MDLAY and MOVE also control the dispatch of stored 
energy. As with the hydro dispatch, if MDLAY is false, then storage units are 
loaded at reduced capacity with conventional units by their marginal cost.* 

If MDLAY is true, then loading is delayed until all the stored energy can be 
discharged at full capacity. The MOVE option controls the storage dispatch 
exactly the same as HYDRO dispatch. 

3. Time- Dependent Units 

As stated previously, time-dependent generation is modeled as a 
reduction in the net load on the utility system. The plant itself is not 
modeled by SYSGEN, so no costs or energy generated are calculated by the 
program. In other words, SYSGEN treats the actual time-dependent plants and 

their output as exogenous inputs to the system. 

The user can simulate increasing penetration levels by having SYSGEN 
multiply any given load reduction curve by integer values up to 99. For 

example, suppose a load reduction curve is the simulated output of a one-MWe 

solar plant, and the user wishes to simulate a system of 5 one-MWe plants. 

The user may specify that five time-dependent units are to be simulated, and 
the program will multiply the input load reduction curve by five. 

SYSGEN is set up to run a base case using the oriainal load curves, and 
then one or more solar cases using the modified load curves. Each set of 

See D. H. Ebbeler, ''Electric Power Generation Systems Probabilistic 
Production Costing and Reliability Analysis," (JPL working paper), April 
1981, for the appropriateness of interpreting this calculated cost as a 
marginal cost and for using that cost in both the hydro and storage 
dispatch algorithms. 
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cases can contain up to four load reduction curves, and there can be a total 
of twenty cases run at one time. The Input data for time-dependent units arc; 

1) Total number of cases to be run. A base and one solar case would 
be a total of two cases, one base and two solar cases would be 
three cases, and so on. 

2) The number of load eduction curves to be Included In each case 
must be specified. The base case would normally use zero load 
reduction curves and a solar case would use from one to four 
reduction curves. 

3) Each load reduction curve corresponds to one type or size of 
time-dependent unit. But as explained above, specifying more than 
one time-dependent unit will cause these curves to be increased 
appropriately. The user must Input the number of units to go with 
each load reduction curve. 

As an example of how this case structure may be used, suppose one wants 
to set up solar plants at two sites and there are three from which to choose. 
We can find the two that are the most complementary by running case sets on 
the three possible combinations of sites. Three load reduction curves, one 
for each site, would be generated by the appropriate solar plant simulation. 
SYSGEN would produce four cases: the base case without the solar output and 

one case for each site combination. Each solar case would include two load 
reduction curves. The best combination of sites would be the case with the 
highest fuel and/or capacity credit. 

The input option MLRED (described on p. 38) is used to specify whether 
time-dependent units will be simulated, MLRED should be set to true if any 
load reduction curves will be used. 

The following is a summary of the input options used to control hydro, 
storage, and time-dependent units. 

4. Hydro and Storage Dispatch Options: MDLAY, MOVE 

MDLAY controls the hydro and storage dispatch strategy. If MDLAY 
is set to false, then reservoir hydro units are always loaded first, at 
reduced capacity to generate all their energy. Storage units are loaded as 
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soon as their marginal costs put them in the loading order.* If MDLAY is set 
to true, then hydro and storage units are delayed until they can generate all 
their energy at full capacity. 

MOVE also controls the hydro and dispatch strategy. If MOVE is set to 
false, limited energy plants are loaded only at valve points of other units. 
That is, tests are made on the viability of bringing up a storage or hydro 
plant only after the previous increment has been completely loaded. If MOVE 
is true, then it is possible to split increments and load the hydro or storage 
unit at points other than the valve points of other units. Setting MOVE to 
true will result in more efficient use of hydro and storage energy, but the 
running time will be longer. If MDLAY is set to false, MOVE is automatically 
set to false. 

5. Storage Cost Option; MSTOR 

MSTOR controls the storage programs. If MSTOR is set to false, 
then the marginal cost of storage is set to the average cost of base loaded 
energy. The expected energy available is taken from the input reservoir 
size. An approximation of base load energy supplied is made on the basis of 
excess base load energy available, disregarding capacity constraints. If 
MSTOR is true, then the storage algorithms are implemented. The dispatch of 
the storage is controlled in either case by MDLAY and MOVE. 

6. Load Reduction Option; MLRED 

MLRED is used to determine if any load reduction curves will be 
used in the simulation. If MLRED is set to false, then no cases with 
time-dependent units will be run. 

Again, see D. H. Ebbeler, ''Electric Power Generation Systems 
Probabilistic Production Costing and Reliability Analysis," (JPL working 
paper), April 1981, for the appropriateness of the non-delay option 
(i.e., setting MDLAY to false). 
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F. LOADING ORDER 

SYSGEN will load units in order of increasing marginal cost, subject to 
the following constraints: 1) valve points of a unit must be loaded in order 

(the third valve point cannot be loaded until the first two have been loaded); 

2) hydro and storage units will be interpolated where they can discharge all 
their energy to minimize costs, subject to two of the input options discussed 
previously, MDLAY and MOVE; and 3) the order may be modified to meet spinning 
reserve requirements, subject to the input option MSPIN. A modification due 
to the spinning reserve requirement is illustrated in the following example. 

Suppose three coal plants called A, B, and C are being loaded. If the 
economic loading order is A, B, C, but the summed capacity of A and B did not 
meet the spinning reserve requirement, and the summed capacity of A and C did, 
then C would be loaded with or before B. The user may specify how far into 
the economic loading order the program may search to find a unit to meet 
spinning reserve (e.g., one unit down the line, two units, all units, etc.). 
There are two loading order options. The first option, LORDOP, is used to 
specify the loading order at the group level (i.e., base, intermediate, or 
peak), and the second, MLORD, is used to specify the entire loading order by 
unit or valve point. 

1 . Loading Option One: LORDOP 

LORDOP is the only input option that is not a true/false input. 

This option will sort plants by marginal cost within groups, where the user 
specifies the order of the groups. For example, one may want all base units 
loaded, followed by intermediate and then peak. Or base and intermediate 
units can be sorted together, followed by peak. This is indicated by a 
three-digit number, XYZ, where: X position = base group 

Y position = intermediate group 
Z position = peak group (1 < X, Y, Z < 3) 
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and the values of X, Y, and Z determine when the group will be loaded. 

For example, 

XYZ » 123: base, intermediate, and peak units are sorted 

separately. All base units are loaded before any 
intermediate units, and all intermediate units are 
loaded before peak units. 

XYZ " 112: base and intermediate units are sorted together and 

loaded before peak units. 

XYZ = 321: base, intermediate and peak units are sorted 

separately. All peak units are loaded first, then 
all intermediate units, then the base units. 

(Default Option: XYZ = 111) 

2. Loading Option Two: MLORD 

MLORD controls the loading order computation. If MLORD is set to 
false, the user determines the entire loading order instead of having it 
computed by the program. Only one loading order may be specified, and it is 
assumed to be the same for all time periods. If a plant is unavailable 
because it is on maintenance, retired, or not yet installed, then it is 
automatically skipped in the loading stack. The capacity of hydro plants and 
storage plants are adjusted so that they discharge as much energy as possible 
at their designated loading point. If MLORD is false and MSPIN is true, 

SYSGEN will compute the cost of keeping the necessary units on spinning 
reserve, but it will not change the loading order. If the units being loaded 
have valve points, the units may be partially loaded. For example* with two 
5-valve point units, the first three valve points of unit A can be loaded, 
then the first two from B, then one from A, and so on. The only constraint is 
that in order to load the jth valve point of a unit, all previous valve 
points must be loaded first. 
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6. OUTPUTS FROM SYS6EN 

SYS6EN has seven report modules, described below. The user can control 
which modules are printed with the five report options explained at the end of 
this section. 


1. Initial Plant and System Report 

This is an echo print of the input data. (See pages A-1 to A40 of 
the sample output file in Appendix A.) 

2. Sorted Limited Energy Plant Report 

The sorted limited energy plant report writes out conventional 


hydro and storage information showing the order in which they are considered 
for loading. 


Variables in report ; 

Hydro/Storage load 
number (ID) 


Unit index 

Hours per time period 


= units' position in the loading stack, 
e.g., if STORAGE ID “ 2, then that 
unit will be loaded only after the 
storage unit with ID “ 1 has been 
loaded, if unit 1 is available. 

= the unit's identifying number. 

= the number of hours that the unit can 
generate at full capacity. (This is 
the reservoir size divided by the 
uni t ' s capac i ty. ) The hydro arrays 
are sorted by this number from 
greatest hours to smallest. 

MWhs per time period = reservoir size of the unit. 

3. Probabi! ity Curve Report 

This report writes out information for the equivalent load demand 
curve. All values (except the area) are in MWs for the curve. (See, for 


example, page A-13 of the sample output file.) 

Variables in report ; 

Load curve spacing = number of MWs between each array 

point on the curve printed beneath. 

Minimum demand = minimum customer demand. 
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Maximum demand 
Equivalent demand area 
Demand curve 


4. Plant Report 

The plant report gives the 
loaded. (See page A-11 in the sample 
Variables in report ; 

Unit index = 

Un it name = 

Unit type = 

Unit valve point = 

MW added = 

Expected start-ups 

Added expected energy = 


maximum equivalent demand. For the 
initial demand curve, this is the 
peak demand. For the final demand 
curve, this is the peak demand plus 
installed capacity. 

area under the probability curve. 

For the initial demand curve, this is 
the energy demand on the system. For 
all other curves, this value does not 
have physical significance. 

equivalent demand probability curve. 
These values are the probability that 
the equivalent demand will be greater 
than or equal to each demand level in 
the load duration curve. The first 
value on the curve is the probability 
that the demand is greater than one 
load curve spacing; the last value is 
the probability that tha equivalent 
demand is greater than the maximum 
original demand. 


information on each unit after it is 
output.) 


input order of the unit (internal 
identification number). 

user identification for the unit. 

class name and loading type, e.g., 
coal, base. 

valve point of the unit currently 
being loaded. 

capacity loaded in this step (MW). 

expected number of times the unit is 
started up. Reported only for the 
first valve point. 

energy (MWh) expected to be generated 
by the current valve point to meet 
customer demand. 




», I 

tf 

k r 
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Fuel cost s 

O&M cost » 

Expected start-up cost * 

Spinning reserve cost = 

Total cost - 

Total capacity factor “ 

Energy used for storage = 

Capacity factor after storage » 

5. 

subperiod after all units have been loaded. 


present worth of the expected fuel 
cost in the tiine period for the 
current valve point (thousand $), 

present worth of the expected O&M 
cost in the time period for the 
current valve point (thousand $). 

present worth of the expected cost of 
starting up the unit. Calculated 
only for the first valve point 
(thousand $). 

present worth of cost of using the 
unit for spinning reserve. Reported 
with the values for the last valve 
point (thousand $). 

sum of fuel, O&M, start-up and 
spinning reserve costs (thousand $). 

unit capacity factor. Capacity 
factor = total energy generated 
divided by the product of the MWs 
loaded and the number of hours in the 
time period. 

energy generated for storage by the 
current valve point (MWh). 

total capacity factor for a base 
loaded plant that is used for storage. 


Subperiod Report 

The system summary report prints out data on the system in each 

(See page A-12 of the sample 


output file.) The first page of the summary report gives the total energy 
generated and total costs for each unit in a format similar to the plant 
report described above. In addition, the total system costs for the subperiod 
and the effective capacity of each unit are printed. The effective capacity 
of a unit is defined as the capacity of a 100%-reliable unit which would 
generate an equivalent amount of expected energy. If the input option MSTOR 


is true, a report on storage losses is written. 
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6. System Report 

The system report prints the following variables for either a 
subperiod or for a time period (see page A-14 of the output); 


Peak demand 

m. 

peak customer demand (MW). 

Customer energy demand 

a: 

original customer energy demand (MWh). 

Load factor 

s 

energy generated/(peak demand x 
hours). 

Unserved energy demand 

a 

expected energy demand which cannot 
be met by the installed capacity 
(MWh). 

Percent energy unserved 

a 

percent of original customer demand 
that cannot be met. 

Loss-of-load probability 


probability that the customer demand 
cannot be met or, equivalently, the 
percent of time customer demand 
cannot be met. 

Magnitude of loss of load 

» 

expected magnitude of each loss of 
load (MW). 

Frequency of loss of load 

a 

number of times in the subperiod or 
period that the load cannot be met. 

Duration of loss of load 

s 

average duration of each loss of load. 

Total expected energy 
generated 

2 

sum of the expected energy generated 
by each unit (MWh). 

Fuel cost 

S 

total system expected fuel cost 
(million $). 

O&M cost 

= 

total system expected O&M cost 
(million $). 

Total cost 

a 

total system expected cost including 
fuel, O&M, start-up, and spinning 
reserve costs (million 5). 

GBtu consumed 

s? 

expected input energy consumed by 
each class, reported only at the end 
of the time period (109 Btu). 


7. Grid File 


This file Is generated by SYS6EN to be used as an Input to SCYLLA. 
SCYLLA computes the worth of tline-dependent units to the system using a static 
economic analysis. The following data are compiled for each subperiod; 

Record Description 

1 load curve spacing 
loss of load frequency 
loss of load duration 

2 loss of load probability 

3 installed capacity for 
each class 

4 energy generated by 
each class 

5 fuel cost for each class 

6 O&M cost for each class 

7 equivalent demand curve (forty 
values centered on a zero loss-of« 
load probability) 

For each time period, the total loss-of-load probability, frequency, and 
duration curves are written at the end of the file. 

The following table shows the options used to control the reports 
printed in the computer output. 

Table 3. SYSGEN Output Options 





Report Option 


Report 

MMAXI 

MAXI MIDI MINI 

M6RID 

1. 

Echo of Input Datfi 

X 



2. 

Hydro/Storage 

X 



3. 

Probability Curve 

X 



4. 

Plant Loading 

X 

X 


5. 

Subperiod Summary 

X 

X X 


6. 

System Summary 

X 

X X X 


7. 

SCYLLA Grid File 



X 
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SECTION IV 
SCYUA 


SCYLLA 1s a linear programming model which performs a static 
optimization of system generating capacity for the solar and no-solar cases. 
Optimal changes In capacity are determined by the equivalent load duration 
curves. Thus, the loss of load probability is held constant In determining 
the solar capacity credit. 

Unfortunately, there are a few problems with this technique. 

Optimization 1s done by picking the cheapest base. Intermediate, and peaking 
class types and then duplicating them to meet the load requirement. With this 
method existing capacity Is not considered, and th? Individual plants In the 
simulation are not modeled. Further, In calculating the value of the solar 
generation, the solar plant should be credited for the capacity that was 
displaced and debited for the new installed capacity. Instead, the capacity 
credit calculated by SCYLLA is the difference between the solar and the 
no -solar optimal capacities. Revision of this methodology Is being studied. 

Most of SCYLLA's inputs are calculated by SYS6EN: class capacities, the 
amount and cost of energy generated by each class, and the equivalent demand 
curve for the system. The only exogenous inputs are capital costs and 
escalation rates by class type. 

Figure 5 illustrates the calculation of the solar capacity credit. In 
the top graph, the y- intercepts are the capital costs in $/MW-year for each 
loading type. The slopes of the cost curves reflect the fuel and n&M costs. 
The bottom graph shows the equivalent load duration curves for the solar and 
no-solar cases. It can be seen that the solar capacity credit for each 
1 oadi ng type 1s B - B • , I - I ' , and P - P ' . 


PRECEDING PAGE BLANK NOT FILMED 


-47- 


0 


100 Percent Time 


P and I are trie original no-solar case optimal capacities of peak and 
intermediate loading units, P' and I' are the nev/ optimal capacities 
for the solar case. Optimal base loading capacity does not change. 


Figure 5. Detemi nation of Solar Capacity Credit 










This questionnaire follows the same qeneral order as Sections II through 
IV of this document and as the interactive input program, PEPS. See Appendix 
B for an explanation of how to set up a data base using batch card input. 

It should be noted that the unit-specific questions, 44 through 61, will 
be repeated for each unit in the simulation. 



I ^ 
.■'U' 



A. TIME PARAMETERS 


1) First year of the study (e.g., 1980); 

2) Last year of the study (e,g., 2000): 

3) Number of time periods In the study (note that this 

should be the last year minus the first year plus V) 

4) Number of subperiods In each period: 

5) Number of "hours" In a time period: 

6) Number of "hours" in each "week"; 

7) Value of the length of the "hour" variable 

(eip., hourly* dally, etc.); 

8) A four (4) character name for the length of the "week" 

variable (e.o., "MNTH" or "WEEK"): 

9) Number of "weeks" In each subperiod: 


Table 4. Weeks per Subperiod 



R. GENERAL AND ECONOMIC INPUTS 

10) Title page heading: A title to appear at the beginning of the 

report which may be up to 120 characters long, including spaces. 


11) Report heading: A title to appear at the top of each page of the 

report which may be up to 40 characters long, including spaces. 


12) In what year's dollars are the input costs? 

(e.g., 1978) 

13) In what year's dollars do you want the output costs 
to be reported? (e.g., 1980) 

14) Conversion factor from input year nominal dollars (12) 
to report year nominal dollars (13): 

15) Consumer price index over the planning period 
(0<CPI<1.0): 


16) Discount rate (real, before tax) to be used in the present 
worth calculation (0<DR<1.0): 


C. ORIGINAL LOAD AND LOAD REDUCTION DATA 


The original load data are the values of electricity demand from which 
the load duration and frequency curves are formed. The load reduction data 
are normally the power output of a time-dependent plant, such as solar 
thermal, photovoltaic, wind or geothermal generation. 

If these data are to be provided on cards, they should be readable on an 
IBM machine and in 12F5.0 format (12 data points per card, starting in column 
21 ). ♦ 

If these data are to be provided on tape, the tape should have the 
f 0 1 1 owi ng char ac ter i sties : 

- standard IBM tape, 9 track 

- 1600 BP I 

- RECFM = fixed block 

- LRECL = 80 

- BLKSIZE = 3200 

- data points in F5.0 format; 12 points per line 

- no label on tape 

Note again that the load data and the load reduction data must contain 
the same number of points, use the same time intervals, start at the same 
point in time, and be in the same units, e.g., MW or kW. 

17) How many data points are in the original load curve? 

(0 < # points < 8,784) 

18) How many load reduction curves are there? _____ 

(0 < # curves ^ 4 ) 

19) How many load duration curves will be calculated?** : 

(0 < # curves < 52) 


Unless FEPS Is being used to create the data base. With the appropriate 
commands, FEPS can transform original data in almost any format into a 
usable SYSGEN data base. For further information, see the FEPS User 
Documentation, JPL Working Paper, November 1980. 

Note that the number of load duration curves is equal to the number of 
subperiods in each period. 
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20) How many load shape points in the load duration curves 

do you want to use? (0 <# points < 100) 

21) FLECTRA automatically creates load duration curves which correspond 
to the subperiod structure of the simulation. Thus, load duration 
curves will be formed sequentially from the original load data, and 
there will be one LDC formed for each subperiod. 

SYS6EN, however, allows the user to assign the LDC's to any 
subperiod desired, For example, a given LDC may be used in more than 
one subperiod, or the LDC assigned to a subperiod may vary hy period. 

In Table 5 on the following page, one should designate in which 
period and subperiod(s) each load duration curve is to be used. Please 
note that if the PEPS program is being used to createthe data base, then 
this information is input in the ELECTRA mode. 


Table 5. Specification of Load Duration Turves 


LOAD 



LOAD 



DURATION 



DURATION 



CURVE 

PERIOD 

SUBPERIOD 

CURVE 

PERIOD 

SUBPERIOD 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 
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D. INPUT OPTIONS 


22) MSPIN = T F (circle one) 

If MSPIN is false, skip questions 23 - 25. 

23) Choose one of the three following methods of stating required 

reserve; 

a) Required reserve is a percent 

of peak load: % 

b) Required reserve is an absolute 

MW value; MW 


c) Required reserve is a fraction of 
the largest unit on line 
(0 < RR < 1) ; ___ 

24) What is the maximum percent of any 
unit to be credited to spinning 
reserve? * % 


25) What is the maximum number of units 

to be displaced in the economic 
loading order in order to meet the 
spinning reserve requirement? 

(integer value) 

26) MULT - T F (circle one) 

(Note; If MULT is set to false, MSPIN is automatically set to 
false.) 


27) 

MFREQ = 

T 

F 

(circle one) 

28) 

MLORD - 

T 

F 

(circle one) 

29) 

LORDOP = 



(3 digit number designating base, intermediate 





peaking order) 

30) 

MDLAY = 

T 

F 

(circle one) 

31) 

MOVE - 

T 

F 

V lire leone) 

32) 

MSTOR = 

T 

F 

(circle one) 

33) 

MLR ED = 

T 

F 

(circle one) 


[ 

I 


E. GENERATION CLASS DATA 


34) Number of generating classes (1<X<34); 

35) Class name of each (can be up to 4 characters long, including blanks), 

followed by class type, where B = base, I = intermediate, and P = peak; 
(e.g., 1. COAL, B . 2. HYPO, P ). 


Table 6. Designation of Generating Class Name and Type 




I 36) In the table below, indicate the capital cost in $/MW of a plant in each 
* class type. Note that these inputs are not necessary if the user is not 
I running SCYLLA. 


Table 7. Capital Costs by Class 





37) How many cost escalation tables will be input? 

(O^ll' tables <10) 

38) How many forced outage multiplier tables will be input? 

(0<# tables <10) 

39) Tables 8, 9, and 10 are all used to assign the cost escalation 
rates and outage multipliers. In Tables 9 and 10, each set of rates and 
multipliers is identified by a table number from 1 to 10. These table 
numbers, or reference numbers, are put into Table 8 to assign the rates 
and multipliers to the appropriate class types. 

Note that there can be only ten total fuel, O&M, and capital cost 
escalation rate tables. 
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Table Reference Numbers for Cost Escalation Rates and Outage 

















TIME 

PERIOD 








F. TRANSMISSION AND DISTRIBUTION LOSSES 


40) For how many Toad levels will time-dependent plant los'ies 
be specified? (1< # levelsiSlO) 

41) In the table below, specify the load levels as fractions 
of the peak load and the marqlnal loss at each level. 


Table 11. Transmission and Olstrlbutlon losses 


G. UNIT SPECIFIC DATA 


If time-dependent units are being modeled, the following data are 


required. 

42) How many cases will be run? (1:2# cases S20) 


43) For each case, specify the total number of load reduction 
curves included, which ELECTRA curves they are (either curve 
1, 2, 3, or 4 if it is a solar case), and the capacity 
multiplier for each curve (i.e., the number of units to be 
simulated.) 


I' 


I 


CASE 


NUMBER OF 
CURVES INCLUDED 


CURVE 

NUMBERS 


NUMBER OF 
PLANTS 


t 


I 




1 


2 


3 


4 


5 


1 

2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
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NUMBER OF CURVE NUMBER OF 

CASE CURVES INCIUDED NUMBERS PLANTS 

6 1 

2 

3 

4 

7 1 

2 

3 

4 

8 1 

2 
3 




NUMBER OF 

CASE CURVES INCLUDED 

W 


15 


16 


17 


18 


19 


20 


CURVE NUMBER OF 

NUMBERS PLANTS 

1 

2 

3 

4 __ 

1 

2 __ 

3 

4 

1 

? 

3 

4 

1 

2 

3 

4 

1 

2 

3 

4 

1 

2 

3 __ 

4 

1 

2 __ 

3 

4 


Unit Index 


44) If this unit is a storage unit; 

a) Charqinq capacity of storage unit (MW): 

b) Forced outage rate for charging cycle 
(0<XS1.0); 

c) “Weekly" energy size of storage unit (MWh); 

d) Generating/charqing efficiency 
(0SX<1.0): 
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} 



Unit Index 

45) If this unit is a conventional hydro plant, then the weekly MWh size of 
the reservoir in each subperiod is required. (Note that these figures 
will be constant over periods, e.g., if the "weeklv" MWh size in 
subperiod 1 is X, then the first subperiod in each year will have MWh 
size X assigned to it.) 

Table 12. Hydro Reservoir Size hy Subperiod 


SUBPERIOn "WEEKLY" MWh SUBPERIOD "WEEKLY" MWh 


t 

► 

I 


k 

I 

I 

I 


! 

I 


! 

i 

I 

f 

i 

IS. 

I 


r 




1 

2 

3 

4 

5 

6 
7 
a 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 


The following Questions must be completed for all conventional units and 
for the hydro, storage, and time-dependent units, as applicable. 

Unit index 

46) Unit name (up to 8 characters)* 

47) Unit class 

(The class name must match one of the class names specified in the 
generation class data, previous section.) 

48) Installment year Week 

49) Retirement year Week 

50) Fuel cost ($/MBtu) 

51) O&M cost ($/MWh) 

52) Cost per start-up ($/start-up) 

53) Cost per MWh to keep unit 
as spinning reserve without 

generating power ($/MWh) 

54) Mean time to repair after 

failure (hours) 

55) Transmission and distribution 

penalty factor (real #^1.0) (default = 1.0) 

56) Total MW capacity of unit 

57) How many loading increments does this unit have? 

(If the unit is either all on or all off, or if 
MULT=F, then the number of loading increments is ''1." 

There may be up to five loading increments.) ________ 

If the unit does not have multiple increments, indicate the average heat 
rate and forced outage rate below: 

58) Average heat rate for unit (MBtu/MWh) _______ 

59) Forced outage rate for unit (0<X<1.0) _______ 

‘ Please note that if the unit names specified are unique, then it will be 
much easier to interpret the computer output. 




¥ 




{■: 

# ¥ 


<S: 






{ • 





% 

ar|i 

I 


3 

1 
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60) Incremental loading data, if applicable: 


Table 13. Incremental Loading Data 




INCREMENTAL 



CAPACITY 

HEAT RATE** 

FORCED 

INCREMENT* 

(MW) 

(MBtu/MWh) 

OUTAGE RATE 



1 


2 


3 

4 

5 


61) Maintenance schedule: 

In the table on the following page, specify one complete maintenance 
cycle by indicating in wivich sufaperiods maintenance is scheduled and the 
number of ''weeks'* per subperiod the unit is unavailable. Please recall that 
the maintenance cycle is referenced to the period the unit was installed, not 
the first period of the simulation. 

Number of years (periods) in the cycle 

(1 < # years S 10) : 


No unit may have more than 5 increments or valve points. 
This is the slope of the total heat rate curve. 
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Table 14, Maintenance Cycle 



5 

[ 


# of Weeks Unit Is Unavailable 


# of Weeks Unit Is 

! 

Unavailable ii 


SUB- 

PERIOD 

PERIOD 

1 2 3 4 5 6 7 8 9 10 

SUB- 

PERIOD 

PERIOD 
1 2 3 4 5 6 

■ 

7 8 9 10 .1 

























! 

I H. LOADING ORDER 


6?) Loading order option HLORD; specifyinq entire loading order (optional): 

Table 15. Loading Order Specification 





1 

i 


UNIT INDEX UNIT NAME VALVE POINT 


t 

i 


LOADING ORDER 

30 

31 







•« '» 


I. REPORT OPTION DATA 

m m 

63) If MGRID is true, then the Grid File Report for a SCYLLA simulation is »• 

vritten. T 

MGRID « T F (circle one) 

Wf IIS 

64) If MINI is true, then the System Summary Report is printed for each 
suhperiod and for each time period, 

MINI = T F (circle one) 

65) If MIDI is true, then all of MINI is printed and the Plant Subperiod and 
Time Period Reports are printed. 

MIDI T F (circle one) ^ 

66) If MAXI is true, then all of MIDI is printed plus the Plant Loading 
Report. 

MAXI = T F (circle one) 

67) If MMAXI is true, then all of MAXI is printed plus the Initial Plant 
Report, the Sorted Hydro Report, and the Probability Curve Report. 

MMAXI = T F (circle one) 

'V V 

i ■ 

r 

■■-Tim 

urn 

mm 

J 

f: 

mm 
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APPENDIX A 

SYSGEN SIMULATION SAMPLE OUTPUT 


The Catalina Island study has been chosen to give a simple illustration 
of the SYSGEN computer model. Peak demand for power on this utility ranges 
from 2.3 to 3,2 MW over the year, with the highest load occurring during the 
summar tourist season. The utility has 5 diesel generators of from 1 MW to a 
little over 1.5 MW capacity. These generators all have the same heat rates, 
forced outage rates, and fuel and O&M costs. 

The following section contains the first subperiod computer output of 
two SYSGEN runs. The first run is the “base case" simulation without solar 
generation; the second run includes a 100 kW solar thermal electric plant in 
the simulation. The energy output of the solar thermal plant was estimated 
with the SES computer model using 1978 insolation data from Barstow, CA. 

The first ten pages of the base case are the echo print of the input 
data. The echo print Is omitted from the second case. Please recall that the 
variables in the computer output are defined in Section III.G. 
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FRACTION OF TIME THAT THE EQUIVALENT DEMAND EXCEEDS THE ARRAY INDEX » SPACING 
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APPENDIX B 
JCL FOR CARO INPUT 


This appendix contains the JCL (Job Control Language) necessary to run 
the ELECTRA, SYS6EN» and SCYLLA simulations on the Caltech IBM computer using 
punched cards. The Input formats of each card set are fully described In the 
original MIT documentations of the three programs; copies of the formats have 
been made and are available upon request. 


A. ELECTRA 


ELECTRA has 5 Input files and 7 output f11es» which contain the 


following data: 



INPUT FILE 

CARD 


Number 

Set 

Input Data 

10 

A.B 

General t economic* and time parameters 

11 

C 

Hourly load data 

12 

D 

Hourly load reduction data 

13 

E 

Transmission and distribution losses 

30 

F 

Generating units data 


OUTPUT FILE 


Number 

Output Data 

22 

ELECTRA debug file 

60 

Base case (no-solar) load duration curve 

61 

Solar case load duration curve 

70 

Base case frequency curve 

71 

Solar case frequency curve 

80 

Base case load-generation correlation 
matrix 

81 

Solar case load-generation correlation 
matrix 


The input card ?ets are punched exactly as described in the input format 
section. The card order for an ELECTRA run is: 

//‘run name* JOB ('account number'), 'programmer name',CLASS>K 
//STEPl EXEC P6M-NEWELCT 
//STEPLIB DO DSN«ACS342.SDDL.L0A0.DISP«SHR 
//FT06F001 DD SYSOUTsB 
//FTlOFOOl DD ★ 
card set A 
card set B 
//FTllFOOl DD * 
card set C 
//FT12F001 DD * 
card set D 
//FT13F001 DD * 
card set E 
//FT30F001 DD * 
card set F 

//FT60F001 DD SYSOUT=B 
//FT70F00: DD SYSOUT=B 
//FT80F001 DD SYSOUT=B 
//FT61F001 DD SYSOUT=B 
//FT71F001 DD SYSOUT=B 
//FT81F001 DD SYS0UT=B 
// 

The 'SYS0UT“B' statement has the output put onto punched cards. If you want 
printed output instead, replace all the *SYS0UT=B* with 'SYS0UT=A*. 
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B. SYS6EN 


SYSGEN has 7 input files, of which 2 are output files from ELECTRA. 

Note that the SYSGEN card sets A, B, and F are similar to the ELECTRA A, B, 
and F sets, but they are not identical. Thus, those inputs canr.jt simply be 
transferred from one program to the other. 

Recall that SYSGEN must be run twice if load reduction curves are being 
used. The first run with the unmodified load curve is called the base case, 
and uses ELECTRA output files 60 and 70. The second run with the modified 
load is called the solar case, and ELECTRA output files 61 and 71 are used as 


the input files. 



INPUT FILE 

CARD 


Number 


Input Data 

10 

A,B 

General, economic, and time parameters 

15 

C/1,2 

SYSGEN load data variables 


C/3 

ELECTRA output file 60 for base case 
run or 61 for solar case run 

20 

D 

ELECTRA output file 70 for base case 
run or 71 for solar case run 

25 

E 

Generation class data 

30 

F 

Generating units data 

35 

G 

Preventative maintenance data 

40 

H 

Loading order specification (input 
only if the entire loading order is 
being specified) 

SYSGEN only has two output files. The content of the first file is 

determined by the report options specified in the intput data (MINI, MIDI, 

etc., in card set A). 

The second 

file contains the grid file for a SCYLLA 

run. If SYSGEN is run 

twice, then 

there will be four output files total: two 


Again, the input cards are punched exactly as described In the SYS6EN 
Input format section. The card order for the two SYS6EN runs Is as follows: 
Base Case (Run 1) 

//•run name' JOB ('account number'), 'programmer name',CLASS«K 
//STEPl EXEC P6M«NEWSY6 
//STEPLIB DD DSN«ACS342.SD0L.L0AD,0ISP»SHR 
//FT06F001 DO SYSOUT-A 
//FTlOFOOl DD * 
card set A 
card set B 
//FT15FOOI DD * 

card set C, using ELECTRA output file 60 
//FT20F001 DD * 

card set D, using ELECTRA output file 70 
//FT25F001 DD ★ 
card set E 
//FT30F001 DD * 
card set F 
//FT35F001 DO * 
card set 6 

//FT40F001 DD * Used only If loading order 

card set H Is specified 

//FT45F001 DD SYS0UT*A,DCB=(RECFM=FBA,LRECL=133,BLKSIZE«133) 

//FT50F001 DD SYS0UT=B 


Solar Case (Run 2) 

//'run name' JOB ('account number'), 'progranmer name',CLASS»K 
//STEPl EXEC PGM-NEMSYG 
//STEPUB DO OSN«ACS34a.SODL,LOAD,DISP<*SHR 
//FT06F001 DD SYSOUT-A 
//FTlOFOOl 00 * 
card set A 
card set B 
//FTlSOOl 00 * 

card set C, using ELECTRA output file 61 
//FT20F001 00 * 

card set 0, using ELECTRA output file 71 
//FT25F001 00 * 
card set E 
//FT30F001 00 * 
card set F 
//FT35F001 00 * 
card set G 

//FT40F001 00 * Used only if loading 

card set H order is specified 

//FT45F001 00 SYSOUT=A.OCB=(RECFMsFBA,LREGL=133 ,BLKSIZE=133) 
//FT50F001 DO SYS0UT=B 
// 
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C. SCYLLA 

SCYLLA's Input files comprise the input card sets A through F for SYSGEN 
(and H if necessary), and the base and solar case grid files output from the 
two SYSGEN runs. 


INPUT FILE 

CARD 


Number 

Set 

Input Data 

10 

A,B 

Identical to SYSGEN 

25 

E 

Identical to SYSGEN 

30 

F 

Identical to SYSGEN 

40 

H 

Identical to SYSGEN 

60 

C 

Identical to SYSGEN 

61 

D 

Identical to SYSGEN 

80 


Base case grid file output from SYSGEN 

81 


Solar case grid file output from 
SYSGEN 


The card order to a SCYLLA run Is: 

//’run name' JOB ('account number'), 'programmer name',CLASS>K 
//STEPl EXEC P6M=NEWSLL 
//STEPLIB DO DSN=ACS34E.SDDL.LOAD,DISP»SHR 
//FT06F001 DD SYS0UT=A 
//FTlOFOOl DD * 
card set A 
card set B 
//FT25F001 DD ♦ 
card set E 
//FT30F001 DD * 
card set F 

//FT40F001 DD * Used only if loading order 

card set H is specified by user 

//FT60F001 DD ♦ 

ELECTRA output file 60 
//FT61F001 DD * 

ELECTRA output file 61 
//FT80F001 DD * 

SYSGEN output file 50, base case grid file 
//FT81F001 DD * 

SYSGEN output file 50, solar case grid file 
//FT45F001 DD SYS0UT=A,DCB=(RECFM=FBA,LRECL=133,BLKSIZE=133) 
//FT50F001 DD SYS0UT=A 
//FT55F001 DD SYS0UT=A 
// 
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In suirmary, the simulation programs have the following Inputs and 


outputs: 

1) ELECTRA Inputs are the card sets A, B, C, D, E, and F specified In the 
ELECTRA Input format. ELECTRA output— files 60, 61, 70, and 71— are in 
the form of punched cards, and comprise the load duration and frequency 
curves for the base and solar cases. 

2) SYSGEN inputs are the card sets A, B, C, 0, E, F, 6. and H specified in 

the SYSGEN Input format. Card sets C and D contain the ELECTRA output 
of punched cards. SYSGEN output Is In two forms: printed output and 

punched cards. The printed output will contain the production costing 
and reliability results of the simulation. The punched cards will be 
the grid file for the SCYLLA program. There will be one printed output 
and one set of punched cards for each SYSGEN run. 

3) SCYLLA inputs are the card sets A, B, C, D, E, and F specified in the 
SYSGEN input format, plus the grid files on the card output from 
SYSGEN. The output of a SCYLLA run is printed. 
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APPENDIX C 

JPL MODIFICATIONS TO SYS6EN CODE 


JPL has identified several problems with the MIT SYS6EN model involving 
both the theoretical foundations of the model and the impleniehtation of the 
theory in the software. Conceptu?! snc- descriptive errors in the MIT 
technical documentation are detailed In a JPL working paper by Dr. D. Ebbeler, 
"Electric Power Generation System Probabilistic Production Costing and 
Reliability Analysis." Dr. Ebbeler's findings are based In part on another 
JPL working paper by Dr. G. Fox, "A Stochastic Formulation of Electric 
Generation System Reliability Measures." 

Listed below are the areas in which Ebbeler and Fox have identified 
problems.* As of this report's publication, not all of these problems have 
been corrected, but the corrections that have been made are noted. A separate 
report is being prepared in which errors in and suggested corrections to 
algorithms and formulas are presented explicitly and illustrated with test 
case results. 

1) The interpolation formula for converting load profiles to load 
duration curves is incorrect. JPL has replaced the algorithm by 
simple linear interpolation. 

2) The deconvolution algorithm is inefficient and unstable for unit 
availability parameters less than 0.5. JPL has replaced this with 
an efficient stable algorithm for those cases. The original 
deconvolution algorithm was also modified to correct for a 
numerical inaccuracy which occurs when deconvolving sufficiently 
small units. 

* Ebbeler, D. H. and G. Fox, "Summary of SYSGEN Changes," Jet Propulsion 

Laboratory, Working paper, March 1981. 


3) The method used by the convolution algorithm to guarantee that the 
sum of expected energies served by the loaded plants and the 
expected unserved energy be fixed at each iteration results In an 
Incorrect computation of all expected energies by distorting the 
peak demand level at each Iteration, JPL has modified this 
algorithm to assure a correct peak demand level at each Iteration, 

4) The algorithm for computing spinning reserve operating costs Is 
Incorrect. JPL has not made corrections, 

5) The algorithm for optimal dispatching of a charged conventional 
reservoir hydro-electric plant was derived In the deterministic 
framework. JPL has not made changes for the probabilistic 
framework. 

6) The storage charging algorithm 1s suboptimal in both the 
deterministic and probabilistic frameworks. The algorithm for 
dispatching a charged storage plant using average cost is 
incorrect. The algorithm for optimal dispatching of a charged 
storage plant was derived in the deterministic framework. JPL has 
not corrected the charging-discharging problem in the probabilistic 
framework . 

7) The treatment of Intermittent generator output as ergodic is 
inappropriate. In the absence of a stochastic model of 
intermittent generator output, JPL models such output as a 
deterministic reduction of customer demand. 


